Optical cooling of atoms in microtraps by time-delayed reflection 
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We present a theoretical analysis of a novel scheme for optical cooling of particles that does not in 
principle require a closed optical transition. A tightly confined laser beam interacting with a trapped 
particle experiences a phase shift, which upon reflection from a mirror or resonant microstructure 
produces a time-delayed optical potential for the particle. This leads to a nonconservative force and 
friction. A quantum model of the system is presented and analyzed in the semiclassical limit. 
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I. INTRODUCTION 

Standard techniques for optical cooling of atoms 
mostly rely on spontaneous emission from the atom to 
carry away the excess momentum. This basic process 
has been demonstrated to be highly efficient for cooling 
of two-level or multi-level atoms and for various laser 
configurations in one, two, or three dimensions, applying 
various polarization states and frequencies 1 -. However, 
a common requirement is for the atom to exhibit a sin- 
gle, albeit possibly degenerate, ground state such that 
no atomic population is lost from the cooling cycle by 
population transfer into internal states decoupled from 
the laser light. These cooling mechanisms are thus not 
well suited for most atomic species and hardly for any 
molecules at all. 

Cavity-mediated cooling mechanisms^ 3 -^^ have been 
suggested to address these shortcomings of laser cooling 
methods. In this case, only a dipole interaction between 
the particles and a near resonant cavity mode is required. 
While this can be fulfilled for a much larger variety of 
particles, resonator alignment and loading of the particles 
into the small mode volume of an appropriate cavity is 
difficult. 

We have recently proposed an alternative method for 
cooling particles that exhibit an electric dipole moment 
but no closed transitions^. In this 'mirror-mediated cool- 
ing' scheme, a laser-driven particle interacts with its own 
image in a mirror and the time-delay incurred during the 
reflection is exploited to introduce a non-conservative ele- 
ment into the dipole force, thereby leading to friction and 
cooling. This time delay may be induced by a delay line, 
e.g., an optical fiber, but more conveniently could arise 
from an integrated optical micro-resonator-Z&^ i 10 ! 11 , pos- 
sibly combined with a plasmonic microstructure for local 
field enhancemen t 12 i 13 . 

Here we investigate the challenges involved in model- 
ing this novel cooling mechanism. We discuss how the re- 
quired processes differ from free-space or cavity-mediated 
laser-cooling methods and propose a theory to overcome 
these new difficulties. We finally discuss the main results 
obtained from our method. 



II. 



CONCEPTUAL DIFFERENCES BETWEEN 
OPTICAL COOLING METHODS 
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FIG. 1: Schematics of optical cooling methods: (a) Free- 
space laser cooling, (b) cavity-mediated cooling of a particle 
between two partially reflective mirrors, (c) mirror-mediated 
cooling of a particle using one perfectly reflecting and one 
partially reflective mirror. 



We start our discussion with a brief review of the fun- 
damental concepts behind three different laser cooling 
schemes: (i) free-space laser cooling, (ii) cavity-mediated 
cooling, and (iii) mirror-mediated cooling. 

Free-space laser cooling methods, such as Doppler 
cooling 14 , polarization gradient cooling 15 , or velocity- 
selected coherent population trapping 1 utilize a number 
of laser beams focused from various directions on a small 
sample of atoms, as illustrated in Fig.[lja). In each case, 
the laser and atomic configuration is set up in such a 
way that the atom-light coupling is velocity-dependent. 
For example, the photon scattering rate or the spatial 
positions where scattering predominantly happens can 
depend on the atomic velocity, leading to an average net 
cooling effect. The excess momentum of hot atoms is 
thereby transferred to spontaneously emitted photons. A 
mathematical description of the cooling mechanism thus 
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contains a number of classical, stationary laser beams, 
and a quantum description of the atomic degrees of free- 
dom. To remove the excess momentum, the excited states 
of the atoms are coupled to a heat bath, which is assumed 
to be memory-less (Markovian)^!. 

In cavity- mediated cooling schemes^^, Fig. [ljb), 
the atoms sit inside an optical resonator and are coupled 
to one or a few resonator modes, while laser beams are 
used to pump the resonator through its mirrors. Dipole 
coupling transfers excess momentum coherently from the 
atoms to the cavity modes, which subsequently decay 
through the partially-transmitting mirrors. Mathemati- 
cally, the classical stationary laser beams are now cou- 
pled to a quantum system comprising the atom degrees 
of freedom and the quantum state of the discrete cavity 
modes. The modes in turn are coupled to a Markovian 
heat bath. 

In mirror-mediated cooling, finally, the particles sit 
outside a coherent delay device, Fig.QJc), which could be 
as simple as a single mirror, a piece of optical fiber with 
an inscribed Bragg grating, or formed by an integrated 
optical resonator on a chipZ&^ i 10 ! 11 . Via the dipole cou- 
pling, the particle imprints a phase on the pump beam 
which returns to the position of the particle after reflec- 
tion with a time delay. Hence, the total dipole potential 
experienced by the particle depends on the interference of 
the incoming beam with the time-delayed reflected beam. 
If the particle has moved during the light round trip, this 
leads to a non-conservative force that can be exploited 
for extracting energy from the particle motion. A math- 
ematical description of this situation therefore requires 
two key ingredients, different from free-space cooling and 
cavity-mediated cooling: (i) the pump beam itself must 
be described as dynamic, and (ii) cooling relies on the 
system state at an earlier time, thus it is non-Markovian. 
On the other hand, no heat bath is required. 

Instead of introducing a system memory in time, we 
may also decompose the field into a continuum of modes. 
The model presented here will follow this latter approach. 
For the sake of simplicity, we model the polarizable par- 
ticle as a single two-level atom and assume that the time 
delay arises from a significant distance (several meters) 
between the atom and the mirror. However, it is envis- 
aged that in a practical realization the time delay will 
arise from reflection by an integrated micro-optical res- 
onator or a similar structure. Moreover, we restrict the 
analysis to a single spatial dimension. The electromag- 
netic modes at angular frequency uj thus are standing 
waves with mode functions f(u,x) = sin(u;a;/c) if the 
mirror is assumed at x = 0. Upon adiabatic elimination 
of the internal degrees of freedom of the atom, and treat- 
ing atomic motion semiclassically, we obtain a continuum 
quantum model governed by the Hamiltonian 

H = J h(u — un)a) (u))a{uu)&u) 

g 2 ff . uoix . U)2X„i, ... . , . 

+n— // sin sin a 1 (LUi)a(uj2 )dcjidcj2 (1) 
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where a(a>) and at(ui) are the mode annihilation and cre- 
ation operators, A is the detuning of the atom from the 
driving laser, and g is the atom-field coupling constant. 
For a two-level atom with transition wavelength A and 
excited state decay rate 2r, this coupling constant is re- 
lated to the mode beam waist w bj*^ 
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where a a — 3\ 2 /(2tt) is the atomic radiative cross sec- 
tion. The force operator describing the action of the field 
on an atom at position x is derived from the Hamiltonian 

as 



AH 



(3) 



We note that the general form of the inter-modal cou- 
pling given by the Hamiltonian (fTJ) is valid for any point- 
like dipole scatterer, and thus will also hold for general 
multi-level atoms and molecules. Moreover, in order to 
calculate the force © we are interested only in the time 
evolution of the field in the vicinity of the particle. Hence, 
it is mainly the relative phase between the mode func- 
tions at this position which is important. We thus con- 
clude that, instead of free-space propagation, any dis- 
persive optical device could be used as a delay element, 
which in particular includes resonant micro- or nanos- 
tructures. 



III. PERTURB ATIVE SOLUTION OF THE 
MODEL 

In the following we apply perturbation theory to the 
model introduced above to derive the basic properties of 
the proposed cooling method. We work in the Heiscn- 
berg picture where the mode operators become time de- 
pendent with the dynamics governed by 
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Next, we apply a semiclassical approximation for the field 
modes assuming that the state of every mode is given 
by a coherent state at all times. This is tantamount to 
replacing the operator a(ui,t) by its expectation value 
a(u>, t). We assume that the atom follows a linear trajec- 
tory x(t) — x + v(t — to) such that at a time to 3> 2r the 
atom is at position x. Here, 2t = 2x/c is the round- trip 
time of the delayed reflection. We then expand the mode 
amplitudes into powers of both the coupling g 2 /A and 
the atom velocity v, 



a(u,t) = a (u),t) 



[ai(w,t) + 1161(0;, t)] 



(5) 
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The zeroth order term in g 2 /A corresponds to the field 
without back action of the atom. It is thus independent 
of v and represents the unperturbed driving laser field, 
which is assumed to be monochromatic, 



ao(u>, t) = A5(lu — ujq). 



(6) 



Here, \A\ 2 gives the pump power in units of photons per 
second. Inserting with ^ into (H|), one can derive 
the first order terms in g 2 / 'A analytically by perturbation 
theory yielding 
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Expressions ©-([2]) can then be inserted into ^ to cal- 
culate the leading terms of the force experienced by the 
atom at position x at time to- Most interesting is the 
first-order term of the force in velocity v, since this gives 
the linear friction force and thus describes heating or 
cooling of the atom in the proposed setup. We find 

o 4 

F v (x, t ) = 2Trhk 2 vT\A\ 2 j^ sm(4k x) (9) 

in leading order of x/X. Here, fco = luo/c is the pump 
wavenumber. Note that F v is independent of time. We 
now define the spatially dependent friction coefficient 
p(x) by the relation F v (x) — —p{x)mv and introduce the 
atomic saturation parameter s — |^4| 2 <? 2 /A 2 at the maxi- 
mum of the standing-wave pump. Together with ^ this 
yields 

it Tik 2 

p(x) = ~4sT— \^rsin(4fc x) (10) 

which has units of s~ x and therefore relates to the inverse 
of the cooling or heating time. 

Fig. [5] shows the dependence of the 1/e cooling time 
for rubidium atoms at a point of maximum friction, 
sin(4fcox) = —1, on the beam waist and for different de- 
lay times ranging from 1 ns to 1 /is. Note that this plot 
assumes constant saturation s and thus the pump laser 
power is assumed to increase linearly with the mode area. 
The figure predicts that cooling times of the order of ms 
can be achieved if the pump is tightly focussed at the 
atom and the delay time is of the order of tens of ns. 
Longer delays lead to faster cooling. These conditions 
are comparable with those explored experimentally in, 
for example, Ref. [lj| 

An important feature of the friction coefficient (fTT)|) is 
its spatial dependence with sin(4fcox). This implies that 
the net friction for an extended spatial distribution is 
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FIG. 2: Cooling time of rubidium atoms at position of max- 
imum friction versus beam waist w for s = 0.1 and r = 1 ns 
(solid line), 10 ns (dashed), 100 ns (dotted), and 1 fis (dash- 
dotted). 
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FIG. 3: Averaged friction coefficient (p(x))$ for a trap cen- 
tered at a point of maximum friction versus trap size 8 for 
s — 0.1, r = 10 ns, and w — 0.5 pm (solid line), w — 1 pm 
(dashes), w — 2 pm (dotted). 



zero. Significant cooling via this method thus requires 
localizing the atom in an additional trap. In practice 
this could be achieved by an additional far-off resonant 
beam either propagating parallel to the mirror or form- 
ing another standing wave superimposed on the driving 
beam for the cooling. Alternatively, on-chip microtraps 
can be utilized in conjunction with integrated time-delay 
reflectors. In the following, we will assume a harmonic 
trapping potential. We can then calculate an averaged 
friction coefficient, defined via the loss of kinetic energy 
over one oscillation, as 



{p{x))& 



-AsT- 
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I sin[4fc a; + 45 sin(T)] cos 2 (T)dT (11) 
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where S is the maximum displacement of the atom from 
the trap center during an oscillation. Fig. [3] shows the 
averaged friction coefficient versus trap size, assuming 
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FIG. 4: Stationary temperature T at position of maximum 
friction versus beam waist w for s = 0.1 and r = 1 ns (solid 
line), 10 ns (dashed), 100 ns (dotted), and 1 /is (dash-dotted). 



that the trap center is at a point of maximum friction. 
As expected, for increasing S the cooling force decreases 
since the atom moves away from the point of optimum 
cooling. For trap sizes S > 0.15 A cooling finally turns 
into heating as the atom spends more time in regions 
of the standing-wave pump where the mirror-mediated 
force accelerates the atom. 

So far, we have neglected momentum diffusion due to 
spontaneous scattering of photons by the atom, which to- 
gether with the friction coefficient determines the station- 
ary temperature achievable with this system. To zeroth 
order in the atom-field coupling, momentum diffusion is 
due to the interaction with the unperturbed pump field. 
We may thus expect diffusion to be identical to that ob- 
served in free-space Doppler cooling, where the diffusion 
constant is given by^ 



D 



(12) 



From this the stationary temperature is obtained as 
D hnw 2 -1 



mpix) t 4cr Q sin(4fcoa;) 



(13) 



IV. SEMICLASSICAL MONTE-CARLO 
SIMULATIONS 

In addition to the perturbative solution of the system 
dynamics, we also performed numerical Monte-Carlo sim- 
ulations. These semiclassical simulations have a number 
of advantages: (i) they allow us to include the harmonic 
dipole trap consistently, (ii) they provide solutions to any 
order in the coupling g 2 /A and in the velocity v, and (iii) 
they include momentum and photon number diffusion. 
On the other hand, the numerical treatment requires us 
to restrict the analysis to a discrete set of equally spaced 
modes of angular frequencies Lo k . 

The corresponding set of equations is derived following 
the approach of Refs. [13 which we only very briefly out- 
line here. Starting from the full quantum master equa- 
tion including quantized atomic motion and a Liouvillc- 
type term for spontaneous atomic decay, a Wigner trans- 
form is applied to obtain a Fokker-Planck equation for 
the joint Wigner function of the complex mode field am- 
plitudes ak and the atomic momentum p and position x. 
In the semiclassical approximation, this Fokker-Planck 
equation is equivalent to the following set of stochastic 
differential equations: 



d.r 
dp 
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i 1 \£(x)4-S* (x) - £* (x)4-£(x)] dt 
da; dx 

-U [£{x)^-£* (x) + £* ( x )^-£(x)] dt 
da; da; 

+k t (x - x )dt + dP, 



(15) 

da fc = iA k a k dt- (iU + ~f)£(x)fZ(x)dt + dA k , (16) 

where £(x) — ^2 k aksin(ujkx/c) is the total field ampli- 
tude at x, xt and fc* are the position of the trap center 
and the trap spring constant, respectively, and 7 and Ua 
are the atomic scattering rate and the optical potential 
per photon, respectively. The terms dP and dA k are cor- 
related stochastic white-noise terms describing the spon- 
taneous redistribution of photons between modes by the 
scattering atom and the subsequent fluctuations in atom 
momentum and modal photon numbers. They are given 
by2° 



where for simplicity we have used the non-averaged value 
of the friction (fTUj). This is a remarkably simple expres- 
sion which, apart from the spatial dependence, only de- 
pends on the delay time and the beam cross section. Fig. 
|4]shows the stationary temperatures corresponding to the 
friction curves of Fig. [2] For realistic parameters these 
simple analytic results predict stationary temperatures 
of the order of mK or even slightly below and thus less 
than an order of magnitude above the Doppler limit of 
141 /iK for Rb atoms. However, we emphasize again that 
the cooling in our scheme is based on the dipole force, in 
contrast to free-space Doppler cooling which relies on the 
radiation pressure force. As such, mirror-mediated cool- 
ing uniquely also works in the far-off resonant regime. 



dP 



dAi 



kv^Jh\£{x)\dW - 



27I— £(x)\dW. 
dx 



(idW + 



(17) 
dW_), (18) 



where dWi (i = 0, +, — ) are independent stochastic Ito 
increments with zero mean, (dVFi) = 0, and unit vari- 
ance, {dWf) = d*21. 

Unfortunately, discretization of the mode frequencies 
also implies that the light field is periodic in time with 
a periodicity given by the inverse spectral mode spacing 
2ir/Aui. It is thus not possible to follow a single sim- 
ulation of the stochastic equations (fT4"]) - (|16p to a quasi- 
stationary state. Instead, we perform averages over short 
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FIG. 5: Momentum of a single-particle trajectory in the ab- 
sence of diffusion, (a) Square of momentum versus time, 
(b) Detail of (a) with fitted turning points (crosses) using 
method 1. (c) Fitting method 2: squared momentum after 
subtraction of a least-square fit with a sine function (solid 
line) and linear fit to the result (dashed). Parameters are for 
Rb atoms, A = -lor, r = 0.25/r, w = 0.7 fxm, Au = 0.1T, 
ujt = 0.5 x 27rF, s — 0.076. The trap is centered at a point of 
maximum friction. 



propagation times for ensembles of different initial tem- 
peratures and derive linear approximations for the rate of 
temperature change dT/dt. Numerically, we found that 
this approach works well, however great care is required 
in the data analysis and fitting routines as outlined in 
the following. 

Figure O shows the principles of our data analysis rou- 
tines for the example of a single particle trajectory, where 
momentum diffusion was neglected for the sake of clarity. 
An atom with an initial momentum of 100 hko oscillates 
in a trap with trap frequency u>t, Fig. OJa). A closer look 
at the maxima of the momentum oscillations, Fig. [5fb), 
reveals the cooling effect due to friction. Note that this 
effect is small, only of the order of 10~ 4 for the period 
of time shown here. It is therefore necessary to remove 
the fundamental oscillation with very high accuracy from 
the simulated data, as it can otherwise easily mask the 
effects of friction. 

One possibility (method 1) is to fit each single oscil- 
lation with a harmonic motion. From these fits the po- 
sitions and momentum amplitudes of the individual os- 
cillations are obtained, as indicated by the crosses in the 
figure. Finally, a linear fit to these data points provides 
an accurate measure of the averaged friction coefficient 
which can be compared to the analytic result (jTTJ) . An 
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FIG. 6: Cooling rate dT/dt versus initial temperature. The 
data points (crosses) and error bars are obtained from Monte- 
Carlo simulations using 10 4 trajectories per initial tempera- 
ture. The solid line is a linear fit to the numerical data. The 
parameters are as in Fig. [5] 

alternative method (method 2) to extract the friction co- 
efficient from a trajectory, or from an ensemble average 
over many trajectories in the presence of noise, is shown 
in Fig. EJc). In this single sine function is fitted 

to the trajectory and subtracted from it. The result is 
a curve where the large amplitude oscillations have been 
removed. A linear fit to this curve reveals the average 
slope and thus the friction coefficient. Both methods, 1 
and 2, work well within certain parameter limits, with 
method 1 in general being more accurate but method 2 
being quicker to evaluate. 

The results of one set of simulations of the full system 
of stochastic differential equations (fTl| - (fTo]) are shown in 
Fig. [6l Simulations were performed in ensembles of 10 4 
independent trajectories, where the initial conditions of 
each ensemble were chosen to represent an atomic cloud 
at a given temperature. Every trajectory was propagated 
for a time 60 T -1 , and ensemble averages of the squared 
momentum were taken as a measure of the ensemble tem- 
perature as a function of time. Finally, method 2 as 
outlined above was employed to extract a linear approx- 
imation to the cooling rate dT/dt. The numerical errors 
related to the finite number of simulations per ensemble 
were estimated by applying the same analysis to sub- 
ensembles of 10 3 trajectories. 

Figure [S] shows that the numerically obtained values 
of dT/dt are well approximated by a linear function of 
initial temperature. This linear behavior is expected for 
Brownian motion with linear friction. For very small ini- 
tial temperatures, momentum diffusion dominates and 
the ensemble temperature increases with time. For large 
initial temperatures, on the other hand, friction domi- 
nates and the ensemble is cooled. The stationary tem- 
perature where the two effects cancel is obtained as 
T = 0.69 ± 0.17 mK from the linear fit to the data with 
Gaussian error propagation. By comparison, the sim- 
ple analytic estimate (fT5|) predicts a temperature of 0.76 
mK, which is in surprisingly good agreement given the 
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number of approximations made in the derivation of the 
perturbative result. 



V. OUTLOOK: APPLICATION TO MICRO- 
AND NANODEVICES 

In the analysis presented above, the delay time was 
provided by a meter-scale distance between the particle 
and the mirror. However, we envisage various routes for 
integration of this cooling scheme into microchip devices, 
as outlined in the following. 

Microresonators. The necessary delay time can 
be conveniently achieved using an integrated optical 
resonato r 7 i 8 i 9 ' 10 i 11 instead of a delay line. For example, a 
10 ns delay requires a resonator quality factor of Q ~ 10 7 , 
which is well within the limits of state-of-the-art micro- 
sphere or microdisk resonators. The friction coefficient p 
will then correspond to the average of Eq. (fTD"]) over the 
delay time r for light exiting the resonator after 1, 2, 3, 
etc. roundtrips. 

Plasmonic field enhancement. In free space, the min- 
imum beam diameter is limited by diffraction to about 
one optical wavelength. However, it has been shown that 
microantennas can vastly enhance local field intensities 
by plasmon effects, and enhancement factors of 300 have 
already been demonstrated^. If the particle could be 
placed inside such a microantenna, the geometric factor 
cr a /(TTw 2 ) would effectively be increased by this factor, 
leading to significantly faster, more efficient cooling. 

Optomechanics. The scheme discussed in this paper 
is concerned with cooling of microscopic particles. How- 



ever, conceptually the proposed cooling method could 
also be applied to larger objects, e.g., micromechanical 
mirrors or oscillator a 22 i 23 ' 24 by appropriate device design 
and scaling of parameters. 

With such generalizations we expect that this cool- 
ing scheme may be exploited in integrated atom chips 2 ^, 
quantum information processors 2 ^, atomic clocks^, or 
interferometric sensors 2 ^. 



VI. CONCLUSIONS 

We have analyzed theoretically an optical cooling 
scheme for neutral particles which, in principle, is not 
restricted to atoms but can be extended to any polariz- 
able species. The scheme exploits the finite delay-time of 
light propagating from the particle to a mirror and back 
to the particle, and is thus non-Markovian in nature, in 
contrast to free-space laser cooling or cavity-mediated 
cooling methods. Finally, we proposed extensions of the 
cooling scheme to chip-based devices, which will open 
the road towards practical applications. 
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